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Abstract 

We show, that SN1987A can serve as an astrophysical laboratory for testing the vi- 
ability of the assertion that a new massive neutral fermion is implied by the KARMEN 
data. We show that a wide range of the parameters characterizing the proposed particle 
is ruled out by the above constraints making this interpretation very unlikely. 

(N : 

1 Introduction 

2 : 

The KARMEN collaboration, which studies neutrino from pion decay at the Rutherford 
laboratory has reported an anomaly in the time structure of their signal, which deviates 
^ . from that expected from the standard decay into muons 0, (2|. One possible interpre- 

\ tation of this, offered in the experimental papers, is the production of a heavy neutral 

fermion with a mass of approximately 33.9 MeV and a specific range of values for the life 
■^j- ■ times. It is the goal of this paper to study the (astrophysical) viability of such a fermion 

CN | using supernova 1987A data. 

' As is well known, the KARMEN detector is located 17.5 m away from the target in 

Q\ ■ the ISIS spallation neutron source. Two pulses of 800 MeV protons, of 100s duration, 

^ ! and 0.32/xs separation, impinge every 20 milliseconds on the target [|J], 0. The time 

Q_i| distribution of ~ 5 x 10 3 events , where energy in excess of 10 MeV has been deposited in 

■ the detector, fits very well the assumption that these are interactions of neutrinos from 
stopped 7r + ( and subsequent /i + ) decays with a moderate, time independent background. 
However, the time bin between 3.1 to 4.1 microseconds has an excess of ~ 100 events over 

■ of the 550 expected from the overall fit.|| 

^ . If this 4.5<t effect is not an accidental fluctuation or an instrumental artifact, its ex- 

planation requires a slow ( (3 ~ ^) relatively monochromatic "messenger" traversing 
the distance L = 17.5m in a time of L/{(3c) ~ 3.6/US. It should then deposit an energy 
E > lOMeV in the detector. Since many (about 10 24 ) neutrons were produced over the 
time, is it possible to explain the anomaly via a neutron messenger? The order of 10 
meter iron shielding corresponds to nla ~ 60 mean free paths of neutron interactions. 
One would then expect the neutrons to be absorbed/ diffused and in any event not to 
come in a well defined time. 

However, the transport of neutrons from target to detector can be dominated by one 
specific crack or external path involving reflections. Even in this case the "neutron = 
messenger" hypothesis may not be viable. Indeed, 

(i) The neutrons emerge from the target with high ( ~ lOOMeV) energies. There is 
no obvious , early moderator which is required to slow the neutrons down to ~ lOOKeV 
energies. Further, there is no obvious mechanism for guaranteeing the required time delay 



L pa th I 1 



1 



with angular brackets indicating average (over many collisions with possible moderation) 
along the crack/external path length, L pa th- 

(ii) The interactions in the anomalous time-bins are distributed uniformly over the 
detector rather than near the front edge or the termination of a putative crack. 

(iii) The neutron shielding around the detector was enhanced in the second phase 
of the experiment. This did not stop the steady buildup of the excess of events in the 
anomalous bin. 

(iv) The strongest argument against the neutron hypothesis comes from the energy 
distribution of the events. Barring the unlikely possibility that the detector had substan- 
tial amounts of 3 He and/or 235 U, the generic (n, 7)reactions of the slow neutrons cannot 
deposit energy in excess of the 8 MeV nuclear binding. If the observed time-anomaly is 
enhanced /weakened by an energy cut E > 15MeV , the case against a neutron messenger 
will be considerably strengthened/weakened. Unfortunately, because of possible excess of 
events in the first energy bin 104- 15 MeV in figure 4 of [|| one cannot make this assertion 
yet. 

For the rest of our discussion we will assume that the KARMEN anomaly has no 
simple explanation. We thus have to address the bold hypothesis ]|, [| that a new neutral 
fermion n° exists. The properties of this, n°, (the "Karmino" implied by the KARMEN 
data) are: 

(i)Its mass is precisely tuned to be 

m n o = m % + - - 5KeV = 33.906 ± 0.005MeV (2) 

We assume that in a (small) fraction (= B t {tt + — > n° + fJ- + )) of the cases, the stopped ir + 
decays into n°+/i + . The monochromatic n , tuned to have (5 = 1/60, will then arrive at 
the KARMEN detector 3.6/is later. 

If some of n° particles decay via n ^ e + e~v or via n°— ► 71/ while traversing the 
detector, an energy > 17MeV (or exactly 17 MeV, respectively) will be deposited. This 
would explain the anomaly if the B r (-n + — > /i + n°) and the decay rate say T n o — ► e + e~u 
satisfy 

(ii) 

Br (ir+ -> n+n°)T(n° e + e~v) = - 7^"° = 2.6 x lO^V 1 (3) 

r(n u — ► e + e v) 

Is such a particle defined by equations (2), and (3) above (and most likely having 
week interactions in matter) consistent with other terrestrial and astrophysical data? 
This issue has been considered before Q. Recent new experimental bounds Q, and new 
astrophysical considerations presented below motivate us to reconsider it. 

2 Terrestrial data and particle physics consider- 
ations 

The particle physics constraints on the properties of the singlet fermion implied by the 
Karmen anomaly (Karmino) have been extensively studied in several recent papers 0. [f| 
|(| . Here we revisit them to set the stage for our ensuing discussion and also to incorporate 
the constraints from the PSI search for the KarminoQ. We also find further constraints on 
the parameters of Karmino within the context of some plausible theoretical assumptions. 
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The striking kinematics of the tt + — ► /i + n° decay at rest, are equally striking in decays 
in flight. Energetic pions yield decay muons moving in the same direction and with the 
same speed. This feature has been used in a recent experiment at PSI-which was inspired 
by the KARMEN anomaly -to obtain the remarkable bound |?J 

B r (7r + -» /x+ + n°) < 2.6 x 10" 8 . (4) 
Equation (3) then implies that 

T-( n o_> e+e -„) < 10 3 s . (5) 

The lack of evidence for sharp ~ 17 MeV energy deposition for events in the anomalous 
bin, and a theoretical bias (the need to use loops rather than tree diagram) suggest that 
the above decay, rather than n° — ► v + 7, dominates. 

It is useful to parameterize the it — ► fi + n° and n° — ► e + e~v processes via effective 
local four-Fermi interactions. 

G^ d (x) l5 V u (x) * n o(x)r%(x) , (6) 

and 

G n0 % (x)T"^ Vd (x) y e (x)r'V e (x) , (7) 

with Lorentz/flavor structure which may be different from those for neutrinos in the 
standard Electroweak model. Still we can compare using just phase space considerations^], 
the expected ratio of the decay rates: T n /T^ (rn n o /m^f G 2 ^/G 2 F Using m n o/m^ « 
1/3, r (l «4x 10 5 s _1 and equation (5) we find that 

G 2 

-p£ > 0.7 x 10~ 6 . (8) 

F 

The local effective Lagrangian implies also the crossed scattering process (at energies 
E v » m 2 n0 /2m e ~ GeV) 



u d + e -»■ e + n° , (9) 

ino appearing in 1 

with crossection 



with Vd the neutrino appearing in the n°decay n° — ► e + e v^, and the process e + e — > rPvd 



a(e + e — > n°zAj) G 2 Q , s 

— - ^ (10) 



a(e + e -> z^) G| 

This last process features in SN1987A and supernovae in general. We therefore need 
to know the constraints on the coupling G n o from laboratory experiments to study the 
atsrophysical viability of the Karmino. 

The specific Lorenz structure of the standard 7r + — > ^l + v (or 7r + — ► n°/i + ) decay only 
mildly affects the decay rate. Hence, 

- r(7r+ r ^ + " V) - *r(* + - »V) < 2.6 x io~ 8 , (ii) 



1 For more detailed discussion that takes into account the general Lorentz structure carefully, see 
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where we included the phase space suppression factor (3 ~ 1/60 for ir + — ► n°fi + . From 
the last equation we find that 

G 2 

-£ < 1.5 x 10~ 6 . (12) 

Very short partial life-time in the 10 _6 s < r(n° — > e + e~fd) < 0.6 x 10~ 3 s range can 
be excluded by the following "theoretical arguments" showing that no new interactions 
stronger than the ordinary "weak" interactions are allowed. Indeed, G n o = Gp^e n a 
implies that 

(\ 5 
e;,, 1 = 0.7 x lO'V^s . (13) 

From this we see that life-times shorter than 10~ 3 seconds would require values of G n o 
larger than Gf- However, in the presently accepted approach to field theories one cannot 
postulate at will new non-renormalizable four-Fermi interactions. The latter should be 
viewed only as the low energy limit of massive boson exchanges. The bilinear vertices of 
the boson exchange tree-diagram, ^B^, which are parts of the interaction Lagrangian in 
a fully renormalizable and acceptable field theory. To date , the only such known bosons 
are the vectorial (spin 1) gluons, photon, and W + , W~ , Z° which generate the strong, 
electromagnetic and weak interactions , respectively. Yet, the exchanges could also be 
scalars generating aS + bP four vertices instead of the (V-A) combination for W^. 

The effective G n o is then of the form gi^/m^z x « ' w ith 5iff2V representing the 
Xe + e~ and Xn Vd vertices and is the mass of the new boson. This should be 
compared with gfy/mfy cx Gf (with the same proportionality constant). Since, gw ~ 0.7 
and we do not wish to entertain g > 1, we need mx < fnw i n order to allow for G > Gf- 
However, we will now have either a charged X~ particle exchange in e~Vd — > X~ — > e~n° 
or (and) a neutral X° in the crossed channel.: e + e~ — > X° — > i>dn°, with mx- or (and) 
m x o smaller than mw- 

The new LEP runs at Wcm ^ 200 GeV strongly exclude new charged particles of 
mass smaller 85GeV. Also, the new X° would have strongly manifested as a new narrow 
resonance in electron-positron scattering in the Tristan and LEP regimes. (We note also 
that for electron-positron scattering at energies larger than the mass of X the latter could 
still manifest if there is an emission of an initial internal bremsshtralung photon.) 

The LEP bounds (N u < 3.01 ) on sequential, I wea k 7^ 0, particles with masses < 40GeV 
imply that n° is a weak-Iso-spin singlet. If the decays 7r + — ► // + n° and/or n° — > e + 
e~Vd do involve also left -handed quarks or leptons then by using SU(2)w rotation the 
processes could be related to their isospin analogues. Thus ir + — > fi + n° would be related 
to 7T° — > vn° which is controlled by the same G^ . The latter may again manifest via 

+ N — > N' + n° in the hot proto neutron star. The latter process is more effective there 
(has a smaller Boltzman factor suppression) than /x + n — ► p + n°. Note, however that 
the charged current process is guaranteed by the very existence of four-Fermi coupling, 
equation (6), with no additional assumptions. 

The 7r + — > /U+n process suggests that, if muonic lepton number is conserved then 
?i°has = 1 and hence the decay neutrino in n ^ e + e~Vd should also have = 1. 
If we make the natural (though not mandatory!) assumption that = then by 
crossing and isospin rotation, we could generate yet another decay mode for the muon. 

H + -> e + v e n° . (14) 
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Apart from a small (~ 40%) phase space correction due to m n o/m^+ ~ 1/3 the rate of 
the new mode is given just like for T(n — > evv) by 



n, + - e W) « §g , (15) 



hence 



r( M + -» e+, e n°) ^ 

- e + w ) ~ G|. _ M ' [ ' 

Since the muon decay rate and the electron spectrum and polarization in this decay have 
been studied carefully over the last four decades, one could deduce that : 

e M < 1(T 3 (if n° -» e+e"^) (17) 
Using equation (8) this implies a lower bound on T r fi^ e + erv : 

r n o — > e + e~u^ > Is (if n° — ► e + e~v^j , (18) 

in addition to the model- Independent upper bound T n o_> e + e - v < 10 3 s implied by the 
PSI upper bound on B r (ir — ► /i + n ). 

To summarize, we note that the PSI experiments in conjunction with some plausible 
theoretical assumptions lead to bounds on the lifetime of n° : 1 s < T n o < 10 3 s. Further- 
more, the strength of the four-Fermi coupling of n° to electrons defined by the parameter 
e n o (see the previous section) has a lower bound e n o > 0.7 x 10~ 6 which will imply lower 
bounds on the production of n°in supernovae environments. We consider the impact of 
this on SN1987A observations in the next section. 



3 Limits Derived from Supernovae, Notably Su- 
pernova 198 7 A 

In the standard scenario for a type-II supernova |J, most (~ 99.7%) of the collapse energy 
is emitted via neutrinos. The latter are trapped in the dense hot proto neutron star for 
few seconds during which they diffuse out. The optical signal is delayed by ~ 1 4- 3 days, 
after the shock emanating from the core traverses the progenitor's envelope. The energy 
of the shock (~ 10 51 erg) manifest as the kinetic energy of the ejected envelope. Only 
about a percent i.e, 10 49 erg manifest in the (bolometric) electromagnetic signal. 

This scenario- largely worked out before the SN87A- explosion has been brilliantly 
verified there. A notable exception is the fact that the progenitor was not a red supergiant 
but a blue supergiant. Accordingly, the optical signal appeared an usually short time of 
~ 3 hours after the neutrino pulse that marked the core collapse. 

Many putative novel, features in particle physics can modify the above scenario and/or 
leave other enduring traces. Hence, considerations of SN1987A data (and supernovae in 
general ) strongly limited [ 10 1 the interactions of axions, right handed neutrinos, (gener- 



ated via m-Dirac an d / or via precession of neutrino spin in magnetic field when the neutrino 
posses an anomalous magnetic moment [11], [12]). Finally electromagnetic neutrino de- 
cays are severely limited. 
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We would like to find the ranges of lifetime r and branching B r allowed by the obser- 
vations from SN1987A. For most part we will only assume that B r /r ~ 2.6 x 10~ 11 s~ 1 
as required by the magnitude of the KARMEN anomaly. In terms of the two four-Fermi 
couplings introduced above this equivalent to 

F 

We note that B r ps 2.6 x 10~ 8 ( the upper bound of the recent PSI experiment) and the 
corresponding r « 10 3 s imply roughly equal e n o , e M ( actually e n o = 0.6 x 10 -6 and e M = 
1.5 x 10~ 6 ). As will be shown below, the existence of the n° s with such properties can 
have quite dramatic effects on supernovae astrophysics. 

3.1 Neutrino emission 

The production of n° in the newly formed hot core can proceed via both the G n o and 
couplings. That is via 

e+e" -> n° + u d , (20) 



with crossection 

and via nuclear scattering 

with crossection 



° « GloEl , (21) 
vN -► n°iV , (22) 



a « . (23 ) 

.0 



The latter process can be viewed as due to tt exchange. At low (< 50 MeV) energies 
this yields an effective local four-Fermi vertex. If the 7r° — ► z/n° vertex cannot be inferred 
from 7r + — > we have instead /z + n — ► n°p with crossection 

a(/i+n - n°p) « GjmJ . (22') 

These processes should be compared with the standard neutrino scattering crossection, 
aw GpE^. In the dense (p > 10 15 grcm -3 , i.e. super-nuclear density) and hot ( 
T > 30MeV) core, the resulting neutrino mean free path is: 

/ = — !— = 100cm. (24) 
no~w 

It implies a total number of neutrino collisions, during diffusion 

N V N collisions ~ f y J > 10 8 , (25) 

where R ~ 30km is adopted for the radius of the very hot core. In a fraction e n fB{m n o /T) 
of the collisions we could have, instead of the usual weak scattering, n° production v a + 
N — ► n° + iV , or p + n — > p + n°. 

Note that the Boltznman factor, fs, reflects the energy distribution of the initial 
particles (neutrinos in the present case). The n° themselves need not be in thermal 
equilibrium. For T ~ 30 MeV the appropriate fg ~ (m/T) 3 ^ 2 e~ m ^l T is ~ 1 and for 
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T~ 10 MeV it is ~ 0.5. When only fi + n — > p + n° contributes, the corresponding 
Boltzman factor is f' B ns 2.7 (E/T) 3/2 e~ m ^ T which for T = lOMeV is ~ 1/4. 

When e n o 3> n° production is dominated by e + e~ — ► n frf. To estimate the rate 
of the latter, we note that all collapse models indicate a sizable electron density in the 
proto-neutron star (^- > 0.2 [|T|]) during the first few seconds. The number of weak 
interactions induced e + e~ collisions per positron during these few seconds is now given 
by: 

/ R \ 2 

Nweak-collisions « h — ) ~ 10 6 , (26) 

where we have used lw = (^eCw)" 1 ~ 1000cm. 

bmce G% ~ ancL ^(e+e--^ g|T ~ e ™ ' n production 

occurs in a fraction of e M /s (or e n ofB) of the 10 8 z^A" (or the 10 6 weak e + e~ ) collisions. 
Altogether, the number of produced n°s is N to tai{n®) ~ N to tai{vi) (10% + 10 6 e n o) > 
20fBN to tai{vi) ~ (20 -r- fyNtotaiivi), In the penultimate step we used x + y > and 
e M e n o = 10~ 12 , and finally substituted /b = 1 (or 0.6) for reactions r±, respectively 
and T = lOMeV. Q 

Note the the above rate of n° particle production (which already tends to exceed 
by 5 -T- 20 the initial number of parent neutrinos) is actually a minimum achieved for 
a particular choice of e^/e„o « A^jy collisions/A^- Weak collisions ~ 100. Generically 
equation (29) above leads to N n o > fN v . with / ~ 10 2 4- 10 6 ! n°particles f\ What 
this apparently paradoxical result, really means is the following. Our assumption that 
neutrinos continuously convert into n° s during their entire ordinary, diffusion time 
tdiff ~ lmeanfreepathN co u/c ~ 0.3s is false. Rather, a sizeable fraction converts into n°on 
a much shorter timescale t convers tuff If <C 0.1s. 

//the n° particles escape on this short timescale, out of the hot proto-neutron star, the 
momentary reaction balance will immediately shift toward the n°s. Catastrophic cooling 
of the core on this short timescale will then ensue. The resulting neutrino pulse would 
be drastically shortened in time and reduced in intensity. The 1MB and Kamiokande 
observation of SN1987A neutrino pulses lasted 5-f- 10 seconds and indicated that the total 
energy emitted via neutrinos was > 3 x 10 53 erg, as expected. Hence, these observations 
by themselves exclude almost the complete range of B r ,r (or e^,e n o) parameters and 
thereby the n° hypothesis itself. 

To complete the argument let us estimate the escape time from the core for any single 
n° particle. The n° particles diffuse out on timescale of t n o di ff « tdiff /100 ~ 3ms, 
once the crossection of the n° s on ambient particles {(7 n o^, cr n o e , cr n o u ) is ten times smaller 

2 Due to the chemical potential in the supernova core there is a suppression of the positron density in the 
supernova core, which is not explicit in our discussion, since the effect cancels 

The suppression factor is ~ (j^J e~^/ T . Taking fi ~ 200 MeV and T ~ 60 MeV, this factor is of order 2 x 10~ 3 . 

There is a more dominant mechanism for tiq production that comes from ls^ + e^ — > e~ + tiq which does not suffer 
from chemical potential suppression. As long as the muon neutrinos are equally abundant in the core, this process 
adds to the first term in N total (n°) ss N tota i{v,) (10% + 10 6 e„o) f B > 20 f B N tota i(vi) « (20 4 5)N tota i{v,), and 
leaves our conclusions unchanged. 

3 If only (22') rather than (20) operates then 10 6 e„o — > 10 6 e n a(f' B / fg) but this decrease the total number of 
the n° particles by merely a factor of 5. 
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than the corresponding o\y (the neutrino or other particles weak crossection) . We will 
therefore make the assumption: 

Al: The nuclear and other crossections of n° particles are (at least ten times) smaller 
than those of neutrinos. 

This assumption is clearly consistent with the value of B r and r(or G 2 and G 2 ) 
required for the KARMEN anomaly. The n° particles are devoid of color. Otherwise 
they would Hadronize and certainly could not penetrate 10 meters of iron. They are 
elecromagnetically neutral and as noted above are also singlets of the full SU(2) x U(l) 
electroweak gauge group. Hence assumption Al is eminently reasonable, though given 
the meager information on n° particles it is not however, mandated in an absolutely 
model-independent manner. 

In all the above discussions, we made the implicit assumption that r n o the actual decay 
life-time of n°is longer than R/c ~ 10 -4 s. Therefore the bounds derived apply only for r 
values exceeding ~ 10 -4 . 

Thus, n°with lifetime in the range 1 s < r n o < 10 3 s would be inconsistent with the 
SN1987A observations. A simple way to see this is to note that the corresonding bounds 
for e n o i.e. 10" 6 < e n o < 1(T 3 are in the range forbidden for general coupling of weakly 
interacting particles to electrons [13]. 

One way to avoid this bound would be to give up the theoretically plausible assumption 
that n° has muonic quantum number and e n o ~ 1 in which case all the produced n°'s get 
trapped and this SN1987A constraint can be avoided. In this case, both the n° and 
have to be sterile neutrinos. In this case also one can derive a constraint on e n o from big 
bang nucleosynthesis. The point is that if e n o is indeed of order one, it will inject sterile 
neutrinos v& into the cosmic soup and they will count as one neutrino species. Thus if 
the BBN constraint on additional neutrino species is much less than one j[(J, then the 
n° lifetime has to be less than 10 -5 sec so that any population of z/^'s injected into the 
cosmic soup by n°decays would be diluted by the QCD phase transition. This implies 
roughly that e n o > 0.5. 

Another possibility is to postulate the existence of exotic interaction of both the n° 
and the v& to quarks so that those interactions trap the n°s in the supernova core. Clearly, 
such interactions must have a very specific property ]l7j i.e. the n° and v& must couple 
only to the isosinglet current so that the 7r + 's do not decay directly into the n° in order to 
be cnsistent with laboratory observations i.e. ir + decay mostly to the /u + no- As shown in 
17], this also helps to avoid BBN constraints. To be a bit more quantitative, we may note 
that if the interaction is denoted by £j n j = < ^j=-n Q Tvd{uTu + dTd), then the n°-nucleon 

scattering cross section is given by a n o N ~ ^(gy + 3g\)G 2 F e 2 E 2 . For an average n°energy 
of 100 MeV, this leads to a mean free path l n o ~ 1.5/e 2 so that trapping in the supernova 
can occur for e > 4 x 10 -3 — 3 x 10 -2 . As has been noted in [17], this value of e is quite 
adequate to suppress the contribution of n°to nucleosynthesis. 

It should however be noted that while this can evade the energy constraints, they 
will be of no help in eliminating the constraint from considerations of shock energy and 
propagation time to be discussed in the last section of the paper. 
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3.2 Galactic e + e population and Early 7 Ray Flash from 
SN1987A 

As indicated above the observed neutrino signal from SN1987A provides a sweeping ar- 
gument against the n° hypothesis. Here we note that the timing, intensity and character 
of the ordinary electromagnetic signals from SN1987A, and from other supernovae, also 
exclude a wide range of r (or equivalently B r ) values. Moreover, commulative effects of 
all past supernovae in the history of the galaxy also lead to strong constraints. The key 
assumption required to enable these constraints is : 

A2: The electromagnetic decay mode n° — ► e + e~v a ' is (one of) the dominant decays of 

n°. 

Note that A2 fails to hold if invisible decays like n° — ► 3u dominate. This assumption 
implies that r(n° — > e + e~Vd) is also (close to) the actual r n o, the lifetime of n°. 

For r > 100s at least 20% of the n° emitted will decay outside the progenitor of 
SN1987A radius (~ 3 x 10 12 cm) where the average velocity of n° was assumed to be 
(3 ~ 0.6. The total number of decay positrons is then N e + ~ N(n° decays for r > R) ~ 
0.3iV n o m 0.2 f B N u . « 4 x 10 54 -=- 2 x 10 56 , where N Vi » 10 57 is the total number of 
emitted neutrinos of any given species (say u n ). If the n° production is dominated by 
/i + + n — ► p + n° rather than by v n + N — ► N' + n° and we use T ~ lOMeV the fs can be 



as small as 10 -3 . Following [19] we note that the decay positrons could manifest in two 
different ways: 

(i) positrons can annihilate on electrons (from decays of different n°particles as well 
as ambient) in the vicinity of the supernova progenitor. 

(ii) Positrons that manage to escape into the interstellar medium, slow down and 
keep accumulating over the galactic history. A steady rate of annihilation yielding the 
monochromatic 512 KeV, 7-Ray line is then established. Dar Goodman and Nussinov in 
a work that preceded the explosion of SN1987A, focused on this aspect. It was found that 
even for fs as small as 10 -3 , an extremely broad range: r v > 10 4 s |^| is excluded by the 
observational limits on the 512KeV annihilation line. The occurrence of SN1987A and the 
lack of any observations of an early 7-Ray pulse, delayed by R/c ~ 100s (R ~ 3 x 19 12 cm 
is the radius of the blue supergiant progenitor of SN198978A), with respect to the neutrino 
signal, strongly limit electromagnetic decay at r ~ 100s. Indeed for r < 100s most of the 
n° s will decay within a distance r ~ R from the progenitor. The resulting electron 
column density n e R w = 10 27 10 30 cm~ 2 , of electrons or positrons suffices 
as n e Ra ann ihu a ti on 1, ensuring almost complete annihilation. This results in a 7 fluence 



^« 3 x (10 5 - 10 8 ) 



5 • in8U™-2 



cm 



corresponding to energy fluence (2-^200)erg cm 2 , for a distance of 50kpc to SN1987A. 
Such a fluence is 6 4- 10 orders of magnitude larger than the observational upper bound ! 
1 

Hence, once N e +/N Ui > 10" 4 i. e. once f B > 10^ 3 and T(n° e + e-v d )/r(n°) rj 0(1) 
(namely A2 holds) we can exclude not only r ~ 100s but even r ~ 10s as well, when the 
number of n° decaying outside the progenitor is suppressed by an additional multiplier 
of e- R l^ « e" 16 ! (when (3 ~ 0.6). 

Next, we turn our attention to even shorter lifetimes and show that also n° decays 
within the progenitor could have dramatic observational imprints on SN1987A. The lack 
of these imprints, severely constraining the allowed r values. 



4 such a large r is required for the particles to traverse red giants envelopes 
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3.3 Shock Energy and Propagation time 



In the standard type-II supernova scenario, the bounce shock from the collapsed neutron 
core propagates outwards and traverses the progenitor mantle and envelop. In the process, 
almost all the shock energy is converted into bulk kinetic energy of the ejected envelop. 
When the shock reaches the surface of the star (actually when it gets to optical depth < 1) 
the heated up stellar gas emits radiation in the UV and optical. This marks the beginning 
of the optical light curve. Thus, for SN19897A, the numerical hydro-dynamical models 
|lq] ) indeed are capable of reproducing the 3 hr time delay between the neutrino 
pulse and the beginning of the optical light curve. The propagation time [l^] can be 
approximated as 

/ R\ ( M ej \ 1 / 2 ( E V 1 / 2 , . 

where M e ,-, the ejected mass, is the mass exterior to the radius where the energy E is 
deposited. 

In the present case, one expects to have in addition to the standard (weakened) bounce 
shock, a shock resulting from the thermalisation (on spot essentially) of annihilated e + e~ 
resulting from the decay of the n° particles to e + e~ that occurs within the mantle or 
envelope. This will occur for lifetime of the n° particles r < 30 s and at a mean radius 
of r = v no T, with a spread in r which is of order r. The energy in this shock is an order 
of magnitude larger than that in the standard bounce shock. The total amount of heat 
energy deposited by the n° s is actually their share in the thermal energy of the proto- 
neutron star core, i.e. E exp ~ 3 x 10 52 erg. This is in contarst with the small fraction of 
the neutrino energy that is converted into shock energy. Thus one expects this shock to 
be both more vigorous and also arrive earlier at the stellar surface, 

Assuming emission at a temperature of 5 -f- lOMeV f3 n o ~ 0.7. Using the initial stellar 



structure given by [18] which is based on the model of [|15|, In the left panel of figure 1 
we present the propagation time computed from equation (27) as function of r. 

We see that for r > 0.03s the propagation time is substantially shorter than the ob- 
servational value of 3hr. In the present framework, a propagation time of 3hr requires the 
progenitor to have a much larger radius which would have rendered it a red supergiant 
while it is known to have been a blue supergiant at the onset of core collapse. Also, the fact 
that the hydrogen-rich envelop would have been pushed out before the (now weakened) 
bounce shock from the core collapse could have reached it, is expected to substantially 
reduce the mixing of hydrogen into the iron rich layers, contrary to the observations. 



3.4 Summary and Conclusions 

We have seen above that the various considerations pertaining to SN1987A, and super- 
novae in general, rule out a large range of possible values of r (B r ). In tyhe right panel of 
figure 1, we present for each set of considerations, the range of r values which are ruled 
out. We also present in the same figure the regions of r which are excluded by the various 
particle physics considerations presented in §2. Some overlap, so that certain ranges of 
t are excluded by more than one argument. The wide range of the ruled out r values 
strongly suggests that the n° hypothesis is not viable. 

We dedicate this paper to our colleague Arnon Dar on the occasion of his 60th birthday. 
We would like to thank F. T. Avignoni, C. Rosenfeld, and S. Mishra for early discussions 
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Shock Energy and Propagation Time 
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Figure 1: Left: Shock propagation time, in hours, as function of the life-time in seconds. Right: 
Ranges of life-times excluded by various effects 

on the KARMEN anomaly and Y. Alster, D. Asheri, A. Kerman, and E. Piasetzky for 
discussions on the neutron hypothesis. The work of all three authors has been supported 
by the Israel-US bi-national fund, grant 94-314. The work of I.G. and S.N. also by the 
Israel national Science Foundation, grant 561/99 and the work R. N. M. is supported by 
the NSF grant no. PHY-9802551. 
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